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ABSTRACT 

Estimates made in the 1970’s indicated that a supernova occurring within 
tens of parsecs of Earth could have significant effects on the ozone layer. Since 
that time improved tools for detailed modeling of atmospheric chemistry have 
been developed to calculate ozone depletion, and advances have been made also 
in theoretical modeling of supernovae and of the resultant gamma ray spectra. 

In addition, one now has better knowledge of the occurrence rate of supernovae 
in the galaxy, and of the spatial distribution of progenitors to core-collapse 
supernovae. We report here the results of two-dimensional atmospheric model 
calculations that take as input the spectral energy distribution of a supernova, 
adopting various distances from Earth and various latitude impact angles. In 
separate simulations we calculate the ozone depletion due to both gamma rays 
and cosmic rays. We find that for the combined ozone depletion from these 
effects roughly to double the “biologically active” UV flux received at the surface 
of the Earth, the supernova must occur at < 8 pc. 

Subject headings: molecular processes; Earth; stars: supernovae: general; supernovae: 


SN1978A; ISM: cosmic rays 
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1. Introduction 

Ruderman (1974) suggested and was the first to study the reduction of stratospheric 
0 3 due to enhanced levels of nitrogen oxides caused by incident radiation from a supernova 
(SN). With rough calculations he found that a nearby (< 17 pc) SN would cause a reduction 
of 0 3 by ~ 80% for > 2 yr from the gamma radiation, and a 40 - 90% reduction in O 3 
lasting hundreds of years from cosmic rays (Ruderman 1974, Laster 1968). Reid et al. 
(1978) reached similar conclusions, while Whitten et al. (1976) found a much smaller 
effect. Schramm k Ellis (1996) present a study similar to Ruderman’s in which simple 
analytical scalings were utilized. Crutzen k Briihl (1996) also examine the problem using 
a time dependent two-dimensional (2D) model, considering only the enhanced gamma ray 
flux. Except for Crutzen k Briihl, these studies used simple theoretical models for the 
SN energy and had, at best, one-dimensional (ID) photochemistry models. Advances in 
computing power and atmospheric models now make a more detailed analysis possible. 
Also, SN1987A observations and recent state-of-the-art calculations provide an estimate for 
the total gamma ray input from a core-collapse SN. 

2. Model 

We model the effects of a core-collapse SN on the Earth’s ozone by inputting the 
expected cosmic and gamma irradiation into an atmospheric model. We consider: (1) the 
relatively short-lived (~ 100 d) gamma rays from the initial blast, and (2) the longer lived 
(~ 20 y) cosmic rays accelerated in the SN blast wave. We examine the effects of SN 
distance and impact angle latitude, and determine the ozone depletion averaged latitudinally 
and globally. We use the NASA/Goddard Space Flight Center (GSFC) two-dimensional 
photochemical transport model (Jackman et al. 1990, 1996) after modifications by Vitt k 
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Jackman (1996) to incorporate energy deposition by solar proton events (SPEs) following 
Armstrong et al. (1989) and galactic cosmic rays (GCRs) from Nicolet (1975). The GSFC 
model is zonally averaged, has 18 latitude bands, and represents altitude with 58 evenly 
spaced logarithmic pressure levels (approximately 2 km vertical grid point separation) from 
the ground up to 116 km, and uses a time step of 1 d. 

The observed gamma ray photon spectrum for SN1987A is 

— = 1.7 x 1CT 3 ( f cm -1 s -1 MeV -1 (1) 

dE VI MeV/ 

(Gehrels et al. 1988) between 0.02 and 2 MeV, lasting 500 d at 55 kpc, for a total energy 
9.0 x 10 46 erg. For ease of modeling we set the incident monoenergetic gamma ray photon 
flux Nf by binning this differential flux into 66 evenly spaced logarithmic intervals from 
0.001 to 10 MeV, for a net energy input of 3.3 x 10 47 erg. For a given distance Dsn we 
scale the empirically observed SN1987A spectrum by (5.5 x 10 /.Dsn) i with Dsn i n P c - 
addition, in our final analysis we rescale our results to a total gamma ray energy of 1.8 x 10 
erg because of the following: SN1987A was unusual in that its progenitor was a blue 
supergiant rather than the more typical red supergiant. A recent three-dimensional (3D) 
smooth particle hydrodynamics (SPH) calculation 3 of a SN with a 15M 0 red supergiant 
progenitor gives - [1.8 ±0.7] x 10 47 erg for models which are initially spherically symmetric. 
In this model, the gamma ray luminosity peaks at t ~ 340 d and is within a factor of 10 of 

3 Aimee Hungerford kindly provided the results of her SN calculations in advance of 
publication. Her initial model, sl5s7b from Weaver & Woosley (1993), was taken 100 s 
after bounce and mapped into a 3D SPH code. To calculate the gamma ray spectra and 
total energies, she mapped time slices from the SN calculation into a 3D Monte Carlo gamma 
ray transport code. 
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the peak for ~ 500 d. Since the energies of the gamma ray photons are so much greater 
than those involved in the atmospheric chemistry reactions, the total input energy is more 
relevant than the detailed SN spectrum and we therefore have done all scalings using 
integrated energy inputs. 

Treating the energy bins as monoenergetic beams of geometric mean energy < E{ > 
and integrating over a range of energies (from i = 1 to 66), the incident photon flux is given 

by 

N° = 8.5 x 10~ 3 [E~°- 2 - E'+i 2 ] cm -2 s _1 (2) 

where the total incident energy flux in the monoenergetic beam at the top of the atmosphere 
is Ff = Nf < Ei >. We attenuate the gamma ray photon flux with altitude via an 
exponential decay law, with the frequency-dependent absorption coefficient taken from a 
look-up table (Plechaty et al. 1981). The beam is propagated vertically through a standard 
atmosphere, which is adjusted later in the photochemistry model for the appropriate 
latitude and time of year. The photon flux remaining in a monoenergetic beam is given 
bv Ni , = N°e ~^' Xj , where or, is the column density (in g cm -2 ) measured from the top of 
the atmosphere and Hi is the mass attenuation or absorption coefficient for the zth energy 
bin (i.e., associated with energy < Ei >). The photon flux, AIV, j (in photons cm s ), 
deposited in the jth layer with energy < Ei >, is the difference between iVty-i and Nij. The 
energy flux deposited in the jth layer by photons of energy < Ei > is Fjj = A Nij < E t > 
(in MeV cm -2 s -1 ). The total ionization rate is the sum over all energies 


9tot, j — 


1 66 F 

1 y' hJ 

35 eV fr{ ’ 


( 3 ) 


where 35 eV is the energy required to produce an ion pair (Porter et al. 1976) and A Z 3 is 
the thickness of the jth slab. The energy deposition versus altitude calculated in this way 
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is not dramatically different from that obtained using a detailed radiative transfer model. 

We estimate the SN cosmic ray inputs from the galactic cosmic ray (GCR) ionization 
rate profiles at various latitudes and solar cycle phase (Vitt & Jackman 1996, Nicolet 
1975). The mean rate of nitrogen atom (N) production at different altitudes and latitudes 
is computed by multiplying the empirically computed rate of ionization (from Nicolet 1975) 
by 1.25 (Jackman et al. 1990). The N production enhances nitrogen oxide amounts as well 
as total odd nitrogen, NO y (e.g., N, NO, NO 2 , NO 3 , N 2 O 5 , HNO 3 , HO 2 NO 2 , CIONO 2 , 
and BrON0 2 ). The local GCR energy density ~ 1 eV cm” 3 is dominated by protons with 
a peak at ~ 0.5 GeV (Webber 1998), implying a flux ~ 7 cm” 2 s" 1 . Each incident cosmic 
ray produces ~ 10 4 * * 7 ionizing secondary particles. The GCR ionization rate inputs (from 
Nicolet 1975) were multiplied by 100 to simulate the charged particle flux from a SN at 10 
pc, and scaled by (10 pc/D S n ) 2 for other distances. For D S n = 10 pc, the locally observed 
GCR flux ~ 0.05 erg cm' 2 s _1 becomes ~ 5 erg cm -2 s~\ yielding a total energy ~ 4 x 10 49 
erg over a 20 y run. The corresponding fluence is ~ 3 x 10 9 erg cm -2 . 

The vertical ionization rate profiles are then mapped onto a 2D circular grid, 
representing a zonally averaged spherical Earth. First a series of ionization rate profiles 
as functions of atmospheric slab number and latitude are calculated for a flat Earth for 
incidence angles of 5°, 15°,..., 85° from the zenith, corresponding to 10° wide latitude bands 
centered on 85°, 75°,..., 5°, respectively, for the illuminated hemisphere. These profiles 
are interpolated to produce vertical average ionization profiles as functions of altitude and 
latitude. The zenith frame daily ionization profiles are mapped into 360 intervals of 1° 

4 We thank David Smith, John Scalo, and Craig Wheeler for kindly sharing results from 

their Monte Carlo atmospheric radiative transfer calculation. The robustness of adopting 

simple energy-dependent attenuation coefficients in terms of the energy deposition versus 

altitude was first shown by Chapman (1931). 
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longitude and 18 intervals of 10° latitude. Zonally averaged daily ionization rate profiles are 
produced by averaging over all longitudes for each 10° latitude band, and then input into 

the GSFC 2D model. 

Three types of calculations of 20 y duration were carried out: (1) a “base” simulation 
with no cosmic or gamma rays, (2) “perturbed” simulations with enhanced cosmic ra> 
levels, and (3) “perturbed” simulations with enhanced gamma rays. The “perturbed” 
simulations were compared to the “base” simulation to assess the atmospheric changes. 
For the cosmic ray trials, the charged particle fluxes are held constant for 20 y, and the 
results are then checked to ensure a steady state is achieved. For the gamma ray trials, the 
irradiation is activated on day 60 (March 1) and maintained at a constant level for 300 d 
(the SN duration), after which it is attenuated abruptly while the model is run for another 
19 y, allowing the perturbed atmosphere time to relax to pre-SN conditions. 


3. Results 


Five gamma-ray simulations were performed for SN impact angle latitude (i.e., latitude 
for which the SN is at zenith) Jsn = —90°, —45°, 0°, +45°, and +90 , taking Dsn — 10 
pc. Three more simulations were carried out for zsn = 0° an d ^sn = 20, 50, and 100 pc. 
We find that after 300 d of simulated gamma radiation input, the calculated changes in 
NO y and O 3 column density depend significantly on *sn- The largest increases in NO y for 
D sn = 10 pc are for z SN = ±90°, and exceed 1800% over the poles. The smallest increases 
occur for i S N = 0°, where the maximum increase is < 900% by the end of the run. Ozone 
depletion accompanying the increase in NO y is greatest for ?sn — i90 , with the zone of 
maximum decrease (~ 60%) at 60-70° S latitude for *sn = — 90 . The smallest decrease 
occurs for z S n = 0°, with the region of maximum decrease (about 40%) at 60-90° S latitude 
(Figures 1 and 2). The maximum annual global average decrease for D S n = 10 pc (Figure 
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3) is 27% in Year 2 for ? S n = 0°; the minimum is 18% for *sn — -90°. Simulations for 10, 
20, 50, and 100 pc show a ~ D ^ trend in ozone depletion, where ~ 1.3 < n <~ 1.9 (Figure 

4) . In all cases, the global average ozone depletion from enhanced gamma rays decreases to 
a few percent by the sixth or seventh year. 

Figure 4 shows an increasing deviation from a D^ law for the ozone depletion as 
Dsn decreases. In other words, although the fluence or net energy input from gamma 
rays varies (by assumption) as D the ozone depletion begins to saturate for small 
D S n- Three processes lead to this deviation: 1) The reaction of N with NO limits the 
production of NO y in the limit where the NO y concentration becomes too high; in effect 
NO y begins to become self-destructive (Crutzen & Briihl 1996). 2) Ozone can be produced 
(rather than destroyed) by enhancements of NO y in the troposphere and in the lowest 
part of the tropical stratosphere (Crutzen & Briihl 1996), a consideration which becomes 
increasingly important for higher rates of energy deposition. 3) Finally, the enhanced NO y 
begins to interfere with other families (e.g., chlorine-, bromine-, and hydrogen-containing 
constituents) that destroy ozone, thus reducing the resultant ozone destruction from those 
families (e.g., this mechanism for production of NO y due to extremely large solar proton 
events was discussed in Jackman et al. 2000). 

For cosmic ray simulations, we find that after achieving the perturbed steady state, 
NO y column densities increase at all latitudes and for all times of the year. For Dsn — 10 
pc, the largest increases in column NO y (650%) occur over the North Pole in late summer, 
while the smallest increases (200%) occur over the South Pole in winter, with an average 
global increase of about 350%. Corresponding column O s levels decrease globally and in all 
seasons by as much as 40% over the North Pole, in late summer and by as little as 5% over 
the equator during May-August for an average global decrease of 22%. Our global ozone 
reductions are less than those of Ruderman (1974) and Reid et al. (1978), but consistent 
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with Whitten et al. (1976). 


4. Discussion 

In order to quantify the impact on life of different levels of ozone depletion, we use the 
results of Madronich et al. (1998). From their Table 1 and Figure 2, we take a factor of ~ 2 
increase in the biologically active UV flux to be a threshold for significant biological effects. 
This increase corresponds to an ozone depletion of ~ 47%, given that the biologically active 
UV flux scales roughly as the reciprocal of the ozone column (Madronich et al. 1998). 
There will probably be an interference between the increases in NO y caused by the gamma 
rays and cosmic rays, therefore a summation of the separately computed ozone depletions 
from the two effects will result in a maximum estimate of the total SN influence. Summing 
our gamma ray and cosmic ray depletions for Dsn = 10 P c > an< l taking into account that 
our adopted energy is larger than that found in the latest SN study mentioned earlier, we 
obtain a fiducial “critical distance” to significantly disrupt ozone of D crit ^ 8 pc for a SN 
with a total gamma ray energy ~ 1.8 x 10 47 erg. 

The impulsive addition of energy into the atmosphere by other types of catastrophes 

5 Although the exponential, or Beer-Lambert model, for attenuation is a good 
approximation at a given wavelength, it is not appropriate for biological effects at the Earth s 
surface. The biological spectral sensitivity depends on the integrated UV flux between about 
290 nm and 330 nm. One can show analytically that the biologically active UV flux varies 
roughly as the reciprocal of the ozone column given that (i) the ozone absorption cross 
section decreases exponentially with wavelength in this spectral region (290-330 nm), and 
(ii) the biological action spectrum decreases exponentially with wavelength over this same 
spectral region (Madronich et al. 1998). 
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may also lead to the production of NO y , therefore the accompanying destruction of ozone 
may not be unique to the SN scenario. A dilution factor of ~ 1(T 22 for the cross section of 
the Earth as seen from 10 pc means that the SN gamma ray energy of 2 x 10 47 erg translates 
to ~ 2 x 10 25 erg intercepted by the Earth. An asteroid with mass ~ 3 x 10 g and speed 
~ 25 km s _1 such as the one responsible for the Cretaceous- Tertiary extinction (Alvarez et 
al. 1980) has ~ 10 30 erg of kinetic energy (KE), so that even if as little as ~ 10 5 of the KE 
were channeled into atmospheric heating, it would rival the SN input. On the other hand, 
the creation of NO y through a sudden event would probably be less long lasting: First, 
the transport of NO y from the stratosphere to the troposphere would occur over a few 
years. Tropospheric NO y exists mainly as HN0 3 , and would be removed by precipitation. 
Second, a localized event in the atmosphere would probably not be communicated globally. 
For instance, if an asteroid generated NO y along its trajectory through the atmosphere, it 
would not necessarily mix globally on a time scale short compared to its destruction time 
scale. A period of sustained energy input is required to maintain continuously high levels of 
NOy if the depletion of ozone is to be effective. 

For a nearby SN to disrupt life on Earth, there must be some finite probability that 
such an event has occurred within the past several hundred million years. For our galaxy 
the SN rate is ~ 1.5 per century (Cappellaro et al. 1999). Clark et al. (1977) estimate the 
Sun should pass within 10 pc of a SN during each spiral arm passage, implying a SN rate of 
r vj 10 Gyr -1 . Two key issues affecting the SN rate are (i) the vertical scale height /isn over 
which supernovae (SNe) are distributed, and (ii) the displacement of the solar system from 
the galactic plane. The study of Clark et al., for example, ignores the vertical disk structure 
and concludes that every passage through a spiral arm results in a close-proximity SN. An 
opposing view is presented by van den Bergh (1994) who adopts /tsn — 300 pc, thereby 
giving ~ 0.3 Gyr" 1 for SNe within 10 pc. Maiz-Apellaniz (2001) utilizes a sample of 0-B5 
stars, precursors to core-collapse SNe, obtained from the Hipparchos catalog to determine 


- 11 - 


h SN = 34 pc. He also finds the Sun to be ~ 24 pc above the galactic plane - a value that is 
fortuitously less than /i SN - The value h SN ~ 30 pc is a factor of ~ 10 less than that used by 
van den Bergh, and leads to a SN rate intermediate between van den Bergh and Clark et 
al., i.e., ~ 3 Gyr” 1 for D SN < 10 pc, or ~ 1.5 Gyr" 1 for D SN < D CTit ~ 8 pc, where we used 
the fact that the rate varies as D | N . 

One potential caveat to this line of reasoning in which the assumed SN fluences vary 
as is the possibility that the Earth may lie within a bubble of hot gas caused by recent 
SNe, so that when a new SN occurs the cosmic ray influence within the bubble may be 
more isotropic and therefore stronger. Benitez et al. (2002) make this argument concerning 
recent SNe in the Local Bubble and the Scorpius-Centaurus OB association which may have 
generated ~ 20 SNe in the last ~ 10 Myr. The cosmic ray flux (f > c r < 1.4 x 10 7 erg cm' 2 
y-i quoted by Benitez et al. using the minimal distance for the Sco-Cen SNe Dsn > 40 pc 
is actually about a factor of 10 less than we adopt in this work for Dsn = 10 PC leading 
us to believe the effect discussed by Benitez et al. may produce ozone depletions of onlv 
a few percent. Nevertheless the complications introduced by having the Earth be inside 
a hot bubble which may enhance the cosmic ray flux during an interval in which several 
successive SNe occur are interesting and worthy of more detailed studies. 


5. Conclusion 

In summary, we have calculated detailed atmospheric models to determine the extent of 
the reduction in ozone due to elevated levels of odd nitrogen induced from gamma rays and 
cosmic rays produced in a local SN. Our procedure is as follows: (i) We utilize a detailed 
2D model for the Earth’s atmosphere, incorporating the latest advances in photochemistry 
and transport, (ii) For the gamma ray spectrum we take as input the observed spectrum 
from S N1 987 A, scaling the total energy to the most recent value determined by workers 
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investigating core-collapse SNe with a red supergiant progenitor, (iii) For the cosmic ray 
spectrum we adopt scaled-up values of the empirically observed ionization rates in the 
atmosphere from galactic cosmic rays, (iv) To estimate the frequency of local SNe we 
utilize recent estimates of global SN rates for spiral galaxies and the results of a recent 
investigation into the vertical spatial extent of core-collapse SN progenitors. Our primary 
finding is that a core-collapse SN would need to be situated approximately 8 pc away to 
produce a combined ozone depletion from both gamma rays and cosmic rays of ~ 47%, 
which would roughly double the globally- averaged, biologically active UV reaching the 
ground. The rate of core-collapse SNe occurring within 8 pc is ~ 1.5 Gyr~ l . As noted 
earlier, our calculated ozone depletion is significantly less than that found by Ruderman 
(1974), and consistent with Whitten et al. (1976). Given the ~ 0.5 Gyr time scale for 
multicellular life on Earth, this extinction mechanism appears to be less important than 
previously thought. 

It is a pleasure to acknowledge stimulating conversations with Frank Asaro, Narciso 
Benitez, Adam Burrows, Enrico Cappellaro, Aimee Hungerford, Sasha Madronich, Frank 
McDonald, John Scalo, David Smith, Floyd Stecker, Sidney van den Bergh, William 
Webber, Craig Wheeler, and Stan Woosley. 
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FIGURE CAPTIONS 

Figure 1. Contour plot of global average of ozone depletion during Year 2 from gamma 
irradiation, for Dsn = 10 pc and 2 sn = 0°. 

Figure 2. Contour plot of global average of ozone depletion during Year 20 due to the 
cosmic irradiation, assuming an enhancement of 100-fold in the empirically observed GCR 
ionization rate, and ?sn = 0°. 

Figure 3. Peak annual average reduction in global ozone from gamma irradiation for 
Dsn = 10 pc as a function of ?sn (average over Year 2). 

Figure 4. Annual average reduction in global ozone from gamma irradiation from a SN 
as a function of Dsni assuming zsn = 0° (average over Year 2). The solid line indicates a 
dependency of 
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